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Abstract — Energy Harvesting (EH) is a novel teclinique to 
prolong the lifetime of the wireless networks such as wireless 
sensor networks or Ad-Hoc networks, by providing an unlimited 
source of energy for their nodes. In this sense, it has emerged 
as a promising technique for Green Communications, recently. 
On the other hand, cooperative communication with the help of 
relay nodes improves the performance of wireless communication 
networks by increasing the system throughput or the reliability 
as well as the range and efficient energy utilization. In order 
to investigate the cooperation in EH nodes, in this paper, we 
consider the problem of optimal power and rate allocation in the 
degraded full-duplex Gaussian relay channel in which source and 
relay can harvest energy from their environments. We consider 
the general stochastic energy arrivals at the source and the relay 
with known EH times and amounts at the transmitters before the 
start of transmission. This problem has a min-max optimization 
form that along with the constraints is not easy to solve. We 
propose a method based on a mathematical theorem proposed 
by Terkelsen 1 1 1 to transform it to a solvable convex optimization 
form. Also, we consider some special cases for the harvesting 
profile of the source and the relay nodes and find their solutions 
efficiently. 

Index Terms — Convex optimization, degraded Gaussian relay 
channel, energy harvesting, , resource allocation. 

L Introduction 

Energy Harvesting (EH) has emerged as a promising solu- 
tion to the perennial energy constraint of wireless networks 
such as wireless sensor networks (WSN) or Ad-hoc networks, 
which have limited battery sources El- Also, EH is devel- 
oped to be used as a foundation of green communication 
networks 13|. This will be more critical as increasing energy 
consumption of highly-demanded ubiquitous networks soon be 
the main cause of global warming. Energy harvesters collect 
ambient energy from the environment or the resources such 
as body heat and convert it into usable electrical energy. 
Conventional devices that harvest energy from the environment 
are solar cells, water mills, wind turbines, microbial fuel 
cells, vibration absorption devices, thermo-electric generators, 
piezoelectric cells, etc. EH nodes have access to an unlimited 
source of energy in contrast to conventional battery-powered 
nodes; however, the Umitation on the EH production rate and 
its sporadic nature necessitates the sophisticated utilization of 
scavenged energy. 

Early works on EH sensor nodes are presented in 
B. In Q, shannon capacity of EH sensor nodes, which 



sends information over an AWGN channel is studied. Also, 
the effects of energy buffer with different condition and also 
inefficiencies in energy storage are studied. Similar study to 
|7| is presented in |8| which showed that AWGN channel 
capacity with stochastic energy arrival (i.e., EH) is the same as 
AWGN channel capacity with average power constraint equals 
to average recharge rate. Also, two achievability schemes are 
presented, namely, save-and-transmit policy and best-effort- 
transmit policy. Optimal packet scheduling problem in wireless 
single user EH communication system, in which energy and 
data packets are stochastically entered the source node is 
considered in ["91. In the aforementioned paper, in order to 
minimize the transmit time of the data packets, transmission 
rate adaptive ly changes according to data and energy traffics. 
The study in f9l is then extended to broadcast channel fTOl, 
1 11], multiple-access channel [12], interference channel [,13]. 
two-hop network lfT4l and fading channels ifTSl . 

Wireless Relay Channel (RC), ever since introduced, is 
used to overcome the challenges of information transmission 
in wireless channels with progressively improving protocols 
proposed by researchers. Cooperative communications based 
on the use of wireless relay nodes is a specific area of research 
in wireless communications that extensively explored in the 
last decade, from many aspects such as information theoretic 
capacity, diversity, outage analysis, cooperative and network 
coding, resource allocation, etc. Also, resource-constraint net- 
works such as WSN can get benefit of cooperation through 
optimal allocation of energy and bandwidth to the nodes based 
on the available channel state information of those nodes (see 
e.g. 1 16 1 and the references therein). 

Some recent studies have considered multi-hop and relay 
networks where their nodes are capable of harvesting ambient 
energy [14^, fTl\, [TSl. In EH two-hop network considered in 
[14], only the relay node can harvest energy, while in two-hop 
networks studied in IITtI both the source (S) and the relay (R) 
nodes are the EH nodes. In (TS^, half-duplex (HD) orthogonal 
RC with decode-and-forward (DF) relay is considered and two 
different delay constraints, namely, one-block decoding delay 
constraint and arbitrary decoding delay constraint up to the 
total transmission blocks, are investigated. 

In this paper, we consider the problem of optimal power 
and rate allocation for EH nodes of a three-node full-duplex 
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Fig. 1. Gaussian relay cliannel with energy liarvesting nodes. 



(FD) degraded Gaussian RC in order to maximize the total 
number of bits that can be dehvered from the source node 
to the destination (D) node in a given deadhne. We consider 
a general model compared to that of the aforementioned 
papers. In our model, there exists a direct link from S to 
D (in contrast to lT4l. ifTTl ) and also we investigate the 
FD mode compared to the HD mode of Li8J. This causes 
a more complicated min-max optimization problem arises in 
our scenario which has not been encountered in the previous 
studies. Our aim is to transform this complicated min-max 
problem to a solvable convex optimization form, using some 
mathematical background. In two steps, we first introduce an 
auxiUary parameter and then use a minimax theorem of [JJ 
to make our problem tractable. Since the online problem that 
assigns rate and power in real-time to the nodes, in our studied 
scenario, is intractable for now, we consider the offline solution 
that the time instants of energy harvesting and the amount of 
harvested energy by S and R are known before the course 
of the transmission. Furthermore, some special cases on the 
harvesting profile of S and R are investigated. The first and 
the second cases have only one node (S or R), which harvests 
energy from its environment, and in the third case, harvesting 
profiles of S and R are the same. These special cases are 
presented to give some intuition of the main problem. 

The remainder of the paper is organized as follows. Section 
introduces the system model, and section Hill formulates the 
throughput maximization problem for the degraded Gaussian 
RC. In section IIVI we provide the optimal solution for the 
degraded Gaussian RC and in section |V] we investigate some 
special cases. Finally, section [VTl concludes the paper. 

II. System Model 

RC models a three-node network, in which the source node 
wants to communicate to the destination node with the help 
of the relay node. 

Definition 1: A (2"^,ri) code for the discrete memoryless 
RC (DM-RC) {Xi X X2,p{y2,yi\x2,xi),yi x 3^3) with four 



finite sets Xi,X2,yi, 3^2 and conditional probability mass 
functions p(j/2 , 2/i 1 2^2 , a^i ) on x 3^2 consists of (i) A message 
set [1 : 2"^] , (ii) An encoder that assigns a codeword (to) 
to each message to e [l : 2"^] , (iii) A R encoder that 
assigns a symbol X2i{yl^^) to each past received sequence 
yl^^ e yl^^ for each i e [1 : n], and (iv) A decoder that 
assigns an estimate rh to each received sequence 2/2 ^ 3^2* oi" 
reports an error The channel is said to be memoryless in the 
sense that given the current transmitted symbols {Xii,X2i), 
the current received symbols {Yu, Y2i) are conditionally inde- 
pendent of the message and the past transmitted and received 
symbols (m, r/-\ Fa*"^). We assume uniform 
distribution of the message over its set. The average probability 
of error is defined as Pe — Pr{M 7^ M}. 

Definition 2: A rate TZ is said to be achievable for the DM- 
RC if there exists a sequence of (2"^, n) codes such that 
lim„^oo -Pi"^ = 0. The capacity C of the DM-RC is the 
supremum of all achievable rates. 

We consider Gaussian RC with EH nodes depicted in Fig. [T] 
The channel outputs correspond to channel input Xi , X2 is as 
follows 

Yi^aX^ + Zu (1) 
Y2 = Xi + bX2 + Z2, (2) 

where a and b are channel gains of S-R and R-D links, 
respectively, assuming normalized channel gain for the S-D 
link, and we have Zi - 7V(0, N), Z2 ^ 7V(0, N). 
Capacity of the degraded RC is as follows |fT9l 

C= max mm{I{Xi,X2;Y2),I{Xi;Yi\X2)}, (3) 

p(xi,X2) 

and the capacity formula for the degraded Gaussian RC 
with power constraint at S, •'^ui'^) — "-^i' ^^'^ 

Ya=i ^2iM ^ "^2, is given by fTll 
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Fig. 2. EH instants and amounts for S and R witli K = 6. 

C(Pi,P2) = min{Ci,C2} = 



Ci = C 



^Piia^Pi - + ^/b^{a^ - l)PiP2 



max {1, a^} Pi 



P2 



N 



otherwise, 
(4) 



which is achieved by Xi - 7V(0, Pi) and Af{0, Pj)- 

III. Problem Formulation 

Our problem is to maximize the number of bits delivered 
by a deadline T from S to D. S and R harvest energy 
at random instants t'~' , ...,t^ and in random amounts 
El,El,...,E^+^ and E^, El, E^+^ , respectively. If at 
some instants only S or R harvests energy, we simply set the 
amounts of the energy harvested by the other one to zero (see 
Fig. |2]i. The interval between two instants S or R or both 
harvest energy from the environment, will be called an epoch. 
The length of i*^ epoch is I' ^ - t'-^ for i = 1, K + 1. 
So, there are a total of -ftT + 1 epoch with = and 
= T — . We consider the offline problem in which 
the arrival times and amounts are known to S and R before 
the start of transmission; therefore, /*, i = 1, iiT + 1 are 
known ahead of time. We find optimal power allocation for S 
and R in order to maximize the rate from S to D with energy 
causality constraints at S and R. This means that energy cannot 
be utilized in S or R before it is harvested in the corresponding 
node. We formulate the problem as follows: 



K+l 

max > ] 

Pl,P2 ^ 



_{Ci(P^,P^),C2(P^,P^)}: 



: max 

Pl,P2 



K+l 

E 



Equation (|6]l denotes that the powers of S and R should be 
nonnegative, (I?) states the energy causality at S and ^ states 
that power consumption at R should not violate its energy 
causality constraint. Finding the solution of the main problem 
is not straightforward as it has the min-max optimization form 
that cannot be separated due to the FD nature of the problem. 
In other word, since R can receive and send the information at 
the same time, in each epoch we do not know which term (in 
epoch i, Ci{Pl, P|) or C2(Pi, P^)) is the minimum and how 
should we assign the powers of S and R to maximize the sum 
of the rates allocated to all of the epochs. Also, observe that 
in (|4|i the condition that specifies the minimum term, depends 
on the optimization parameters of the maximization problem, 
i.e.. Pi and P2. 

IV. Optimal Solution for Degraded Gaussian RC 

In this section, we propose a method to make the problem 
tractable and try to solve it. We can rewrite the problem by 
introducing < A < 1 as follows 



K+l 



max min \ XC, (P^ 



s.t. 



k 



PI > 0, P2* > 0, 

^pir<Y,E\, 

i=l i=0 
k fc-1 



'l,p;) + {i~x^)c,{Pi,p^)} 

(9) 

1, (10) 



i=l 



i=0 



i = l,...,K- 
k = 1,...,K 



1, 



(11) 
(12) 



where. 



if 
if 

arbitrary if 



CiiPlPi)>C,{Pl,Pi), 
Ci{PlPi)<C^{PlPi), 

Ci{Pi,pr) = c,{Pi,p^). 



(13) 



We then use the following corollary to change the order of 
min and max operators which is the application of a min-max 
theorem of Terkelsen |[T], presented in ll20l . 

Corollary 1 ( [20! Corollary 1]): Let A^; be the d- 
dimensional simplex, i.e. > and X^iLi = 1 ■ 
Let V be the set of probability distributions p{u). Let 
Ti{p{u)), i — l,...,d be a set of functions such that the set 
A, 



CiPi Pi) -4={(ai,a2,...,a<j) e R"* : < Ti(p(u)) for some p{u) e P}, 

(14) 

(5) is a convex set. Then, 



s.t. Pi > 0, P2* > 0, 



fc-i 



J2pII'<J2eI 

i=l i=0 
k k-1 

Y^Pir<Y.Ei, 



i^l,...,K + l, 
k = l,...,K + l, 

k = l,...,K + 1. 



i=l 



i=0 



(6) 

sup min KTi{p{u)) = min sup y ^ \iTi{p{u)). 
(1\ p{u)ev ^^^d ^^^d p{u)ev ^^-^ 

^ ' (15) 
Now, we consider our optimization problem defined in (|9). 
(-g) Theorem 1: Optimal power and rate allocation for the 
source and the relay node in the degraded Gaussian RC 



i=l 
K 



V 



fe-1 



K / k 



fe-1 



. /max|l,a2|pf 
+ (1-A')C( ^ ' 

K+l K+l 

+ E '^^^i + E ^^^2, 



(30) 



fe=l \i=l 



i=0 



fe=l 



problem 

K+l 

min max V (A'd(P|,P*) 
{V} {P-}.{P^} ^ 



(l-y)C2(P^P^)} 



with energy harvesting nodes, is the solution of the following C']^(P]^, Pj)^ T2{p{u)) — C'2(A'^2) ^i^^ c? = 2. Note that 

as we fix the input distribution of S and R to be Gaussian, 
i.e., Xi ^ A/'(0, Pi) and X2 ^ 7V(0,P2), we can replace 
Ti{p{xi, X2)) with Ti{Pi,P2) in the above corollary (we 
rewrite p{xi,X2) = g{Pi,P2))- Now, assume that (01,02) € 
(16) A and (61,62) e A. It means that ai < C^{P^,P2),a2 < 
C0(Pi,P2) as well as 61 < C?(Pi,P2), 62 < Cl{P„P2). 
Then we choose (ci,C2) ~ (jyai + (1 — rj)bi,rja2 + (1 — 
i)b2), < 77 < 1. It is clear that ci < C5'(Pi,P2) 
and C2 < C2(Pi,P2). Hence, we have (ci,c2) G ^. This 
completes the proof. 



s.t. 



PI > 0, P2* > 0, 

k fc-1 

J2pIi^<J2ei 

fe fc-1 

E^2^'<E^2, 

i=l i=0 

< A* < 1, 



i = l,...,if- 



1, (17) 
-1, (18) 



:1,...,K+1, (19) 

1 K 



Now, we decompose our problem into the following two 
1 (20) problems. 

Corollary 2: The problem defined in (ISTT i can be decom- 
Proo/.- We show that the corollary [T] is appHcable to our po^ed into the following problems: 



problem. In other words, we have 

K+l 



K+l f _ 

'V^' • D^^ , /I \^^/=./™D^^\ (Problem 1) :/* ({A'}) = max E A'C'i(Pi%P|) 

^^max^^ min |A*C'i(Pi%P|) + (l-A*)C'2(Pi%P2*)| = ^ {Pl},{Pi} 

' +(1-V)C'2(P^P2^)} 



min max . f A'(7i(Pi\ P|) + (1 - A0C2(Pi\ P^)]. 

{A-}{P^'},{P|} ,=1 L J 

(21) 

Now, we show that every convex combination of the points 
inside A is also in A. Therefore, ^ is a convex set. Consider 
the following mutual information terms 

/(Xi; ri|X2, Q) = H{Yi\X2,Q) - H{Yi\Xi,X2, Q) 
(fc) 

< H{Yi\X2)-H{Yi\Xi,X2,Q) 

= H{Yi\X2)-H{Yi\Xi,X2) 

= I{X,;Y,\X2), (22) 



(a) 



s.t. 



PI > 0, P2* > 0, 

fc fe-i 

E^^r<E^^^ 



1 K + l, 



(24) 
(25) 



+ (26) 



i=l 
k 



4=0 

fc-1 



E^2^'<E^2, fc = l,...,i^+l. (27) 



i=l 



i=0 



(Problem 2) : min /*({A*}) 



s.t. 



0<X'<1, 



1, 



(28) 
(29) 



I{Xi,X2;Y2\Q)'^' H{Y2\Q) - H{Y2\Xi,X2,Q) 

H{Y2\X,,X2,Q) 
H{Y2\Xi,X2) 



< H{Y2) 



'^H{Y2) 



id) 



I{Xi,X2;Y2), 



(23) 



which can be solved, separately. Problem 1 is a convex 
optimization problem as its objective function is concave and 
its constraints are affine, and can be solved by efficient convex 
optimization methods to find its unique maximizer Problem 2 
is a combinatorial problem that can be solved efficiently, too. 
For the first problem, we can write the Lagrangian function 
where, (a) and (d) follow from the definition of the mutual for any > 0, /Xfc > 0, t?fc > and Tyfc > as ^ in top 
information, (b) follows from the fact that conditioning does of this page together with following complementary slackness 
not increase the entropy, and (c) follows from the fact that Q conditions 
is a function of Xi and X2. This completes the proof. 

Also, we can specially show that the aforementioned corol- 
lary is applicable to the Gaussian case. We set Ti{p{u)) = 



fc-1 



a E^^''"E^i =0' k = l,...,K, 



(31) 



fc-1 



i=0 / 

N 

J2^,Pl=0, i = l,...,K + l, (33) 

i=l 

AT 

Y^Tj^P^^O, t = l,...,K + l. (34) 



N 



This problem can be solved by taking the derivatives of the 
Lagrangian function with respect to Pi and P2 and setting 
them to zero, and doing some mathematical manipulation; 
however, the closed form expression for Pi and P2 give 
not any explicit idea about the optimal power assignment 
algorithm. Hence, in the next we prove some lemmas about 
the properties of the optimal solution. We use these lemmas 
in the next section to find the optimal solution for a special 
case. 

Lemma 1: In an optimal policy, transmit rates and powers 
of the S and R are constant within an energy harvesting epoch 
and only potentially change at energy harvesting instants. 

Proof: As we know, Ci ~ gi{Pi,P2), C2 = 92{Pi) 
are nonnegative, strictly concave and monotonically increasing 
function of their variables Pi and P2- We prove this lemma 
by contradiction. Assume that there is a i* S such 
that S and R use A and P2 in {p-^,t*) and A and A in 
{t*,P), respectively. Hence we have 



({t--t^-')Pi+{t'^ 



P - 1 



—gi{Pi,P2)< 



p - p 



(35) 



Similarly 



t* - P- 

P - P- 



-92{Pi) 



P - t* 



92 



p - p- 



T32(Pi) < 



(36) 



Therefore by equalizing the transmitted power within an epoch 
we can reach to a higher throughput. Hence, changing the 
transmitted power of S and R within an epoch is suboptimal. 

■ 

Lemma 2: Whenever the power of source or relay changes, 
it should only increase. 

Proof: This is also due to the concavity of Ci = 
9i{Pi,P2), C2 ~ 92{Pi) and the fact that postponing the 
transmission of energy or shifting it to the right (in energy 
consumption diagram) does not violate the energy causality 
constraint and on the other hand, due to the concavity of the 
rate function in terms of power, more bits per joule can be sent 
by setting the power to a constant value. Therefore, the power 
of S or R never decreases in time, i.e., P^ < A^ — A^ — ' ' ' 
and A <Pi<P2< - -- ■ 

Corollary 3: In the optimal policy, if power of S or R 
changes in an instant, the total harvested energy in the previous 
epoch of that node has been consumed completely by this 
instant. 



V. Some Special Cases 

In this section we present some special cases, which are 
interesting from the practical viewpoint. We can consider cases 
in which only one node, i.e. S or R, can harvest energy 
from the environment. Also there may be scenarios that the 
harvesting process of S and R are the same. These special 
cases do not have the complexity of the main problem and 
presented here to give intuition to the main problem. Although 
three cases are solvable, we only state the solution for the third 
one, which is general compared to the others. 

A. Only Relay Harvests Energy 

In this scenario, the topology of the network is such that 
only R can harvest energy from its environment and S has 
solely a non-replenishable battery. This is equivalent to the 
case that in our system model we set A 7^ 0, E'^ = ... = 
El = 0. In this case we have PI ^ ^ ^ P,yi and the 
capacity formula is as follows 



CliP2) 



c 



(^Va^P-b-^P2+b^/(a^-l}P2) ^ 



c 



if P2 < (a2 - 1)P 
o.w. 
(37) 

As we can see in (|37] |. the condition that specifies the 
capacity formula is only the function of P2 and therefore the 
complexity of the main problem is not exists here; so, finding 
the solution of this problem is straightforward. 

B. Only Source Harvests Energy 

This scenario is in contrast to that of the previous one, 
in which only S can harvest ambient energy while R has 
a conventional non-rechargeable battery. This means that in 
our model we set A ^ 0,i;| ... = = 0; So, 

P2 = ^ = P, V« and the capacity in this scenario is as 



CiiiPi) = <^ 



C 



c 



(Va^Pi-b^P+b^(a2-l)p) ^ 



if Pi > 



N 



o.w. 

(38) 



The same conclusion can be made as the previous case. 

C. Same Harvesting Process for Source and Relay 

This is also an interesting case that both of S and R 
can harvest the ambient energy. In this scenario, S and R 
are considered in the vicinity of each other; therefore, we 
can assume the same harvesting profile for them, i.e., the 
harvesting instants are the same and the harvested amounts 
are scaled version of each other Mathematically speaking, we 



have E] 



2' 



1, 



,K, and thus P^ = ■yPl, \fi, for 



some positive constant 7. The capacity formula is given as 

(^a2-b27+6^7(a2-l))'Pi \ 



Chi (Pi) = 



C 



C 



> 1 or a > \/7+T 



o.w. 



(39) 



We obtain, the Lagrangian of our problem in this case as 

K+l 

mPl}. {P2}, A^, ^, ^)=E min CliPlP^)} I 



VI. Conclusion 

In this paper, we investigated the optimal power and rate 
allocation for a three-node full-duplex degraded Gaussian relay 
channel with energy harvesting source and relay nodes. The 
original problem has a complicated min-max form that is 
not easy to solve. We transformed it to a tractable convex 
optimization problem, which can be solved efficiently. Also, 
some special cases on the harvesting profile of the source and 
the relay nodes were considered. 




(40) [5] 



By setting the derivative of the Lagrangian, with respect to S 
power sequence, to zero, we get the optimal power sequence 
of S as 
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1 




Ki 


1 


1 


2A, 


K2 



PI 



where A, = £ ^ - f 

k—i 

K2 = max{l,a^}/N. 

As we have assumed that El 
> and P^ > 0. 
conclude that P{ > 0, 



a > V7+T, 
a < V7+ 1, 



(41) 



and 



Pi' 



> and E2 > 0, therefore 
Also, according to Lemma 2, we 
Vi e {1, ...,K + l}. This together with 
complementary slackness means that we have i?* = 0, Vi G 
{1, ...,K + 1}, and hence Ai = X^sLi Cfc- Note that according 
to dH]) we have < Pi < ■ ■ ■ and hence P2 < Pi < ■ ■ ■ ■ 
This proves Lamma 2 for this problem. Now, we present the 
following lemma for this problem. 

Lemma 3: In the optimal policy, the power at S or R 
changes only when their corresponding energy causality con- 
straints are active. 

Proof: Here we give a proof for this lemma based on 
complementary slackness condition. Observe that if PI ^ 
Pl'^^ then, according to (l4Tl l. we should have 7^ 0. Then, 
according to the complementary slackness condition of (ISTT l. 
we should have X]j=i Pi^'' ^ Sj=o P'l ~ ^- T^^^^ completes 
the proof of lemma [3] for this problem. ■ 
Now, we can present the optimal solution form as follows: 



arg 



mm 

_i<i<A; 



pv 



A^j=o„_i 1 



P2=i Pi- 



rn 



(43) 



Once the powers of S and R are determined, their correspond 
ing rates can be obtained using capacity formula 
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